Study of the triangular lattice tV model near x — 1/3 
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We study extended Hubbard model on a triangular lattice near doping x = 1/3, which may be 
relevant for the recently discovered superconductor Na x Co02 • J/H2O. By generalizing this model to 
TV fermionic species, we formulate a meanfield description in the limit of large N. In meanfield, we 
find two possible phases: a renormalized Fermi liquid and a a/3 x a/3 charge density wave state. 
The transition between the two phases is driven by increasing the nearest neighbor repulsion and is 
found to be first order for doping x = 1/3, but occurs close to the point of the local instability of 
the uniform liquid. We also study fluctuations about the uniform meanfield state in a systematic 
1/N expansion, focusing on the residual interaction of quasiparticles and possible superconducting 
instabilities due to this interaction. Upon moving towards the CDW instability, the increasing 
charge fluctuations favor a particular /-wave triplet state. (This state was recently discussed by 
Tanaka et ai, cond-mat/0311266 ). We also report a direct Gutzwiller wavefunction study of the 
spin-1/2 model. 
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I. INTRODUCTION 

Motivated by the recently discovered 
superconductivity^ in Na^CoC^ • J/H2O (x rs 1/3) 
and the unusual electronic properties^ of the Na^OoC^ 
series, we began a study of he tV model on a triangular 
lattice^ focusing primarily on the renormalized Fermi 
liquid regime. We pointed out that nearest neighbor 
repulsion can lead to significant renormalization of the 
effective hopping amplitude, which seems to be the case 
in these materials. We also argued that strong repulsion 
drives this model into a a/3 x \/3 charge ordered state 
near commensurate dopings x — 1/3 and 2/3. 

We continue this study here and perform a system- 
atic slave boson meanfield analysis of the Fermi liquid 
close to the charge density wave (CDW) state for dopings 
near x — 1/3. We also consider the residual interaction 
between the quasiparticles and study which supercon- 
ductivity channels are favored in the tV model on the 
triangular lattice. Of particular interest here is possible 
enhancement of some channels due to charge fluctuations 
upon approaching the CDW phase. 

Such studies of superconductivity due to residual in- 
teraction in the models with strong local repulsion are 
familiar in the high-T c field. Scalapino et al. A studied the 
Hubbard model on a 3D cubic lattice in a random phase 
approximation (RPA) and found that d-wave pairing be- 
comes attractive close to the spin density wave transi- 
tion. Kotliar and Liu£ studied the infmite-C/ Hubbard 
model on a 2D square lattice in a systematic large-A^ 
treatment and found that the residual interaction from 
the no-double-occupancy constraint favors d-wave super- 
conductivity close to half filling. More recently, McKen- 
zie et alfi applied the analysis of Kotliar and Liu to ex- 
tended Hubbard model on the square lattice at quarter 
filling and found a transition to the v2x v2 CDW, which 
is driven by the nearest neighbor repulsion. Merino and 
McKenziei studied superconducting instabilities of the 
Fermi liquid near this transition. 

It is of interest to perform similar studies on the 



triangular lattice. RPA treatment in the spirit of 
Scalapino et alX has been done by Tanaka et al& very 
recently. They find that in the regime of interest for 
Na^CoC^ • 1/H2O, a particular /-wave triplet channel is 
favored close to the CDW instability. 

Here, we report a Kotliar-Liu type analysis for the tri- 
angular lattice, which has not been done so far. The 
advantage of our approach over the RPA is that the treat- 
ment of the strong on-site repulsion is better controlled. 

In our study, we also find that the preferred supercon- 
ducting channel near the transition to the a/3 x a/3 CDW 
is the /-wave triplet state discussed by Tanaka et al&. 
This state has the lobes of large A oriented towards M 
points of the Brillouin zone (BZ) boundary as depicted 
in Fig. ^ The origin of this result can be understood 
from the following Fermi surface nesting argument. 

In an RPA- type treatment, the effective interaction at 
wavevector q is 



V cS (q) = 



V(q) 



(1) 



1 + 2 X0 (q)V(q) 

where Xo(q) is positive. This can be also written as 

V eS (q) = V(q) - 2 X (q)V(q) 2 (2) 

where x(<z) = Xo(s)/[l + ^XoivW (q)] is the full suscepti- 
bility in RPA. Now, V(q) is negative over some portion of 
the BZ, in particular at the a/3 x a/3 ordering wavevec- 
tors such as Q — 47r/(3a) x, so for sufficiently large V 
there is an instability at which x(Q) diverges. Eq. (0) 
shows that close to the instability the effective interac- 
tion becomes strongly attractive at the ordering wavevec- 
tors. If there are sections of the Fermi surface that are 
nearly connected by such wavevectors, then we expect 
enhanced superconductivity in channels that utilize this 
nesting. Our calculation using slave boson theory also 
produces attractive effective interactions due to charge 
fluctuations. The difference is that now the effective po- 
tential depends on k and k' of the electrons, rather than 
q = k — k'. Near the Fermi surface we write this as 
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FIG. 1: Schematics of the nesting argument for doping x = 
1/3. External hexagon shows triangular lattice Brillouin zone, 
while internal hexagon shows reduced BZ for the v3 x \f?> or- 
der. Flat surfaces of the reduced BZ boundary are connected 
by the \/3 x \/3 ordering wavevectors perpendicular to the 
surfaces. Near these momenta transfer, V c a(9,0') becomes 
attractive close to the CDW instability. We also show the 
/-wave triplet gap function that we find to be the dominant 
superconducting channel from the residual interaction; this 
channel utilizes the near nesting of 9m points indicated with 
black diamonds. 



model with N fermion species. Near x = 1/3, two can- 
didate states are studied — a uniform Fermi liquid and a 
state with V3 x %/3 charge order. The phase diagram is 
established. We then study fluctuations over the uniform 
saddle point in 1 /N expansion, and focus on the residual 
quasiparticle interactions from such fluctuations. This 
approach can be viewed as a version of RPA that treats 
the on-site constraint systematically, and provides a more 
quantitative justification of the above nesting argument. 

We conclude this Introduction with one remark. Ear- 
lier worksi2iii*i2iiii4 considered t J model on the trian- 
gular lattice and found that the dominant superconduc- 
tivity from the J interaction is d + id singlet state. This 
conclusion was also reached in our earlier report in which 
we studied tJV model and treated J as the main residual 
interaction. The tV model with no J terms predicts /- 
wave triplet state near the CDW order, and offers a way 
to distinguish between these two different pairing mecha- 
nisms. The experimental situation regarding the pairing 
symmetry of Na^OoC^ • J/H2O remains inconclusive, with 
controversy still surrounding whether there is a jump in 
the Knight shift below T c i^i& 



II. EXTENDED HUBBARD MODEL ON 
TRIANGULAR LATTICE 



V c h(8,8') and we find that V e g has attractive compo- 
nents when 9 and 9' are connected by an ordering vector 
Q. 

Figure^shows the situation for doping x = 1/3 (drawn 
to scale) . It turns out that at this filling the Fermi surface 
lies fairly close to the reduced BZ corresponding to the 
\/3 x order. We also show the gap function of the 
dominant /-wave state. This has large absolute values in 
the directions of the M points of the full BZ boundary 
and signs as shown — positive sign near black diamonds 
and negative sign for white diamonds. With an attractive 
Vcs(Q), pairing between the black diamonds and also 
pairing between the white diamonds are very favorable, 
since the black diamonds are nearly connected among 
themselves by ordering wavevectors, and so are the white 
diamonds. 

We should of course consider the near nesting of other 
Fermi surface segments as well. For example, take the 
point 9k and the point opposite to it on the Fermi surface 
shown by the filled circles in Fig.^ These points are also 
nearly connected by an ordering vector and this pairing 
would favor a spin singlet state. As an example, an s- 
wave state gains from all points 9, 9' that are connected 
by V e g(9,9') < 0. The reason why the s-wave channel 
is still disfavored is that there remains significant overall 
repulsion over non-nested V c s(9,9') that it cannot avoid. 
The above /-wave state appears to be best for the overall 
V cS . 

The paper is organized as follows. We first perform 
a slave boson meanfield treatment of a generalized tV 



We consider extended Hubbard model in the limit of 
large onsite repulsion, i.e., the following tV Hamiltonian 

Ht v = - Pg *v c L c j<y p G + ^ V v Ui n i ( 3 ) 

ij ij 

with nearest-neighbor repulsion Vy = V. Pg projects out 
double occupation of sites. The band is less than half- 
filled, with the average fermion density of 1 — x per site, 
and we specifically consider the case t > 0. See Refs.0y 
for a more detailed discussion of the possible application 
of this Hamiltonian to the Na x Co02 system. Here, we 
are primarily interested in the doping near x = 1/3. 



A. Meanfield formalism 

Our meanfield treatment follows closely McKen- 
zie et al£>. The general formalism is the same except 
for an arithmetic difference in Eq. 1141 and some minor 
differences in the analysis. We apply this formalism to 
the triangular lattice tV model. Slave boson formulation 
is used to treat the no-double-occupancy constraint. We 
write c ic — f^bi, and the slave boson Hamiltonian acts 
in the Hilbert space with j\ a Uo + b\bi = 1; the boson 
field b keeps track of the empty sites. 

To formulate meanfield description and also in order 
to go beyond the meanfield in a systematic manner, we 
consider a generalized model with TV fermionic species. 
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The slave boson Hamiltonian is written as 



H = - 



NS 



ij 



which now acts in the space 



fJ i „ + blb i = NS 



(4) 



(5) 



(our large N limit is thus like a thermodynamic limit in 
the fermion flavors). At the end of the calculation, we 
will put N = 2, 5=1/2, and our specific choices when 
defining Hfb are such that this will reproduce the slave 
boson Hamiltonian for the spin-1/2 model. 

Proceeding as in McKenzie et al& and Kotliar and 
Lit£, we write the path integral in the radial gauge 17 



(the spin index a runs from 1 to N). For convenience, we 
introduced parameter S, which is kept fixed as we take 
N — > oo. We also used a particular form for the repulsion 
term. To study the system behavior with doping, we fix 
the total fermion number to NS(1 — x) fermions per site 



Z = / UfDfDrDX exp 



/ driL(r) 
J o 



(6) 



The imaginary time Lagrangian 



i ij ij i i 



To obtain this, auxiliary field Ai(r) was introduced to 
enforce the constraint; the phase of the original boson 
field bi was gauged away, while the absolute value is now 
a real field fj = \bi\ [the corresponding measure is T>r = 
Yli T r<(r)dri(r)]. Chemical potential /i sets the correct 
fermion density. 

We now formally integrate out the fermions and seek 
saddle points of the resulting action in terms of the fields 
ri and Xi. The natural scale for the r fields is r 2 ~ N, and 
upon such rescaling (which we do not perform explicitly) 
the saddle point analysis can be cast into a formal large 
N procedure. 

The saddle point conditions read 

(fiafur)+^ = NS, (8) 

i.e., the constraints are satisfied on average, and 
^ r . 

= iXl = 2NS ^ f- ^(fwfc^ + *i*(/jo-^")] 

3 1 
3 

We seek time-independent saddle points. However, the 
fields may be spatially varying to allow for possible charge 
inhomogeneity. In this case, tPi+zWs £j Vij(NS — r|) is 
an effective potential on site i (cf. Eq.Q), while -wgfirjtij 
is an effective hopping amplitude in such meanficld. 
Among all saddle points, we are to take the one that 
minimizes the free energy specified by Eq. Q. 

We first consider the uniform saddle point, which has 
r 2 = b 2 = NSx and 

*Ai = V = Efe + \v )m) ■ (10) 

k 



Here, tu = Ylr' trr'e~ lki * r ~ r ' and similarly for Vk (we 
will often drop tildes when the meaning is unambiguous); 
L d is the number of lattice sites; /(£) = l/(er* + 1) is 
the Fermi distribution; and ^ is the quasiparticle energy 
measured relative to the Fermi level 

= -xh + ¥ + ^o(l - x) - n . (11) 

The chemical potential is tuned so that 

^E/^) = ^ 1 -- t ) • ( 12 ) 

k 

The uniform saddle point represents a renormalized 
Fermi liquid with effective hopping xtij. We recognize 
such hopping energy renormalization as coming from the 
configurational constraints imposed by the no-double- 
occupancy condition. At this level, the repulsive interac- 
tion leads only to the shift in the bottom of the band. 



B. Fluctuations over the uniform state 

Away from half-filling and for small V, the uniform 
saddle point has the lowest free energy. We can estab- 
lish the region of its local stability by considering small 
fluct uatio ns above the uniform state. Proceeding as in 
Refs. l r lGl we obtain the following quadratic action for 
the fluctuations iXi = tp + idXi, r\ = b(l + Sri) 
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(13) 



where w n is bosonic Matsubara frequency, and the "inverse RPA propagator" is given by 

b 4 n ^ m +q ) ~ /(^) 



2b 2 N 

T rr (q,w n ) = 2b 2 <p- —j-j^2{t k 

k 



(TVS) 2 L d ^ iw n + & - 



r rA (g,o; n ) = r Ar (g,o; n )=i 
L\\(q,U n ) = Tri. / >~ TT c 



26 2 b 2 N ^ /(g fe+g ) - /(&) 



(*k + ifc+ 9 + V q ) 



(15) 
(16) 



The above expressions coincide with Eq. (17) in Ref. 6 
upon replacements t/(NS) -> t/N, V/(NS) -> 2VyA, 
except that the Vo term in T rr replaces Vjt there. 
We can write more compactly 



2b 2 



r, , = ^- V - JxgyiWXo + 2^*1 + X 2 ) , (17) 



NS 



(V g X +X 1 ) 



NS 

r rX = r Ar = i 2b 2 

Taa = X i 
where we used Eq. I|10|) and introduced 



(18) 
(19) 



X p (q,u n ) 



N 



, \p /(£&+<?) /(£fc) / on \ 

tk+ q r- — — r — t — [20) 



(21) 



As discussed in Ref. 0, the local stability of the saddle 
point is determined by the condition det T = T rr T\\ — 
^r>X\r > 0, which translates to 



v ,'4& 2 2Y 
Xo 1 NS Vq + NS 



NS NS 



Xi 
NS 



> . 

(22) 

We now focus on the stability to static perturbations, 
setting bj n — 0. Since Xq > and V q is negative in 
some portion of the Brillouin zone, there is clearly an 
instability for sufficiently large V. 

The local stability analysis specialized to N = 2, S = 
1/2 (and in the zero-temperature limit) for all dopings is 
summarized in Fig. [21 

For x < 0.14 there is an instability at q — for small 
V, which persists as long as 

detr(g = 0) = const x [1 + 2t F u{e F ) + V u(e F )] < , 

(23) 

where we used X o (q=0) = v(e F ), etc., v(ep) is the spin- 
ful density of states per site at the Fermi level. The 



instability is towards phase separation into hole-rich and 
hole-poor regions. This can be seen by examining the 
energy of the uniform state (measured per site) 



H(x) = E t {x) + ^%{\~xf 



(24) 



where E t (x) is the hopping energy per site in the free 
fermion problem with hopping amplitude tx and fermion 
density (1 — x). After some analysis, the condition i|23|) 
is seen to be equivalent to H"{x) < 0, which indeed 
leads to phase separation. Note that moderate nearest- 
neighbor repulsion V stabilizes the uniform saddle point 
against phase separation, and the corresponding bound- 
ary is shown in the lower left hand corner of Fig. 

The instability of the uniform state for large V is found 
to always occur at the V3 x V3 ordering wavevector. 
The corresponding critical (V/t) c is shown with the dot- 
ted line in Fig El In the region designated renormalized 
Fermi liquid the uniform state is stable towards static 
fluctuations at any wavevector. 



C. ^xV3 CDW saddle points 

However, local stability does not guarantee global sta- 
bility of the saddle point. As a specific example, we con- 
sider x = 1/3. The above local stability analysis gives 



V c (x= 1/3) = 2.11 . (25) 

The uniform saddle point has the energy per site 

£wform(z - 1/3) - -6tx{f}j ja ) + 3V(1 - x) 2 (26) 

with (fjpfja) = 0.337 for the uniform hopping at this 
doping. 

On the other hand, consider a competing state with 
complete -\/3 x order which has all charges on the B 
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FIG. 2: Meanfield phase diagram of the tV Hamiltonian. 
For x < 0.14, the uniform state is unstable at q — (to- 
wards phase separation) for small repulsion V, but is stabi- 
lized for moderate V . For large V, the instability occurs at 
the Vo x Vo ordering wavevector and the critical (V/t) c is 
shown with the dotted line. In the region between the dot- 
ted lines, the uniform state is locally stable at all wavevectors. 
The circles show the actual meanfield transitions with general 
v3 x \/3 ansatz; filled circles indicate first-order phase tran- 
sitions, while open circles indicate second-order transitions. 



isfy b 2 A + b 2 B + b 2 c = 3a; . Similarly, we allow differ- 
ent ip a an d <Pb — '■Pc- F° r convenience, we can take 
<PA + Pb + Pc = 0, since there is also the chemical po- 
tential degree of freedom which is tuned to obtain the 
correct fermion density for each set bAiPA- Now, for 
each bA we tune ipA until b\ + (f Acr fAcr) = 1 is satisfied, 
and finally adjust b A so that the other self-consistency 
condition is satisfied. We find all such saddle points and 
compare their free energies obtained from Eq. 10. The 
result is shown in Fig. [21 with circles. Filled circles indi- 
cate the situation similar to x = 1/3, when the transition 
is first order. For open circles, we tentatively conclude 
that the transition point coincides with that determined 
by the local stability analysis and the transition is second 
order. 

We considered in detail the uniform and the -s/3 x V3 
CDW saddle points. Since we are interested in the doping 
regime near x = 1/3, where we do not expect some other 
state to enter the competition, the presented meanfield 
analysis is complete. 



III. RESIDUAL SUPERCONDUCTING 
INSTABILITIES 



and C sublattices, bs = be = 0, while the sublattice A 
is not occupied, bA = 1- The energy per site is 

*W3(z = V3) = K (27) 

It is simple to see that the latter becomes energetically 
preferred over the uniform state for V = 2.00, i.e., be- 
low the local instability point of the uniform state, which 
means that the transition is first order. The following de- 
tailed analysis finds that the the actual transition occurs 
at v; lstordor (a; = 1/3) = 1.98 to a state which is close to 
the above state with complete CDW order. 

To study such possibilities in detail, we consider sad- 
dle points that have the \/3 x V3 ordering pattern. The 
procedure is as follows: We select the A-sublattice of 
the three sublattices, and seek saddle points that have 
bi = b A on the A sublattice and bi = bs = be on 
the B and C sublattices. Of course, these must sat- 



We now aim to go beyond the meanfield description. 
The free energy of the uniform state at next order ml/N 
can be obtained from the quadratic action S^ 2 ' , Eq. i|13[l . 
In terms of the fermionic quasiparticles, this contains 
further effective mass renormahzation as well as resid- 
ual interaction of quasiparticles mediated by the bosons. 
We focus on the residual interaction, in order to decide 
which pairing channel is favored at low temperatures for 
our triangular lattice with hard repulsion. 

The vertices coupling the fermions and the bosons can 
be obtained by examining the Lagrangian Eq. Q , while 
the boson propagators are given by the inverse of the ma- 
trix T(q,u) n ) in the quadratic action, e.g., D rr (q,u n ) = 
(5r(—q,—u) n )5r(q,w n )) = T\\/ detT, etc. From the ex- 
pressions for T and using b 2 ~ N, all propagators are 
0(1/N). 

For our study, we single out interaction terms relevant 
for the BCS instability, obtaining 



b 4 

V eS {h,k 2 ) = D xx - T^gyDrAtk, + t k2 + V q \ 2 + 



b 2 



where all boson propagators are at wavevector q — tonian 
k% — &2- The interaction is of order 0(1/N). If we ignore 

the frequency dependence of the boson propagators, the ~ 1 rr rr, r. \ A A f f 

relevant terms correspond to the standard pairing Hamil- ^pairing - l^jj. VcS t* 1 ' k 2)h u aJ-k u a>J-k2^'Jk 2 ,v 

fci ,fc 2 



(29) 
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FIG. 3: Residual quasiparticle interaction relevant for the 
pairing instability (cf. Eq. I29H plotted for the momenta fci 
and on the Fermi surface. One angle is kept fixed while 
the other is varied over the FS as indicated in the insets, which 
also show the triangular lattice BZ and the reduced BZ for the 
V3xV3 CDW order (see Fig. [TJ. The data is for x = 1/3, 
and the uniform state becomes unstable at V c = 2.11. The 
vertical scale is in units of t. Top panel: The fixed angle is 
in the direction of the K point of the full BZ boundary. As 
we approach the CDW instability, 180° scattering amplitude 
becomes attractive, which is associated with the fact that the 
corresponding Fermi surface points (black circles) are nearly 
connected by one of the y/3 x \/3 ordering wavevectors. Bot- 
tom panel: The fixed angle is in the direction of the M point; 
the three nearly nested points occur at 120° from each other. 



From now on, we fix N = 2, S = 1/2. We first consider 
in detail doping x = 1/3. We are primarily interested in 
the scattering processes when momenta k\ and &2 lie near 
the Fermi surface. Therefore, we visualize the effective 
interaction by fixing the momenta on the Fermi surface 
and plotting V^h (6*i , 6*2) as a function of polar angles. 
This is shown in Fig. |3 where we see how V c ff evolves for 
increasing V approaching the CDW transition. Observe 
that there is residual repulsion even when V = 0, which is 
entirely due to the no-double-occupancy constraint, and 
that this repulsion already has some momentum space 
features. As we increase V from zero, the initial effect is 



to add more repulsion at zero momentum transfer. In- 
creasing V further and approaching the -\/3 x y/3 CDW 
instability, there is a dramatic development in the fea- 
tures associated with scattering at the ^/i x y/3 ordering 
wavevectors. As discussed in the introduction, the quasi- 
particle interaction at such momentum transfer becomes 
attractive close to the critical point, and this is clearly 
seen in the plots. In Fig. [3 the critical point V c = 2.11 is 
defined from the quadratic fluctuation analysis. As dis- 
cussed earlier, in the meanfield, we find instead lst-order 
phase transition at a somewhat lower V^ lstordcr = 1.98. 
However the features in V e s(0, 9') are already enhanced 
at the actual transition since it occurs close to the insta- 
bility point and has significant critical fluctuations, so we 
will mostly ignore the distinction in what follows. 

We now study BCS instabilities due to this residual 
interaction; the analysis below is valid for either singlet 
or S z = triplet superconducting channels. Wc follow 
Refs. HI and define coupling constant associated with 
each channel A k ~ g a (k) 



da k 
\v k \ 



da k 



(2ir) d \v kl 



Tg a (k)*V eS (k,k')g a (k') 



da k 



\g a (k)f 



where the integration is over the Fermi surface elements 
da, and v(k) = Vfc£(fc) is the Fermi velocity. For an 
attractive channel, we must have c a > 0, and the tran- 
sitions temperature is given by the weak coupling BCS 
expression (which is appropriate in the present case) 



T c [a] = f2exp[-l/ccj 



(31) 



the frequency cutoff is roughly O ~ xt, since the energy 
integration is over the entire band. 

Representative triangular lattice tight binding har- 
monics that cover main symmetry classes are listed in 
Table [i] Four of the ansatze have their real-space A rr i 
nonzero on nearest-neighbor bonds only, the NNN /-wave 
ansatz has next-nearest-neighbor bonds, and the i-wave 
ansatz requires even further neighbor bonds. Since V e fi 
is real, only real ansatze need to be considered. In- 
deed, in this case the problem of finding a harmonic 
with the largest coupling constant is equivalent to a real 
symmetric eigenvalue problem; the easiest way to see 
this is to discretize the da integrals and define ip a {k) = 
g a {k){da k /\v k \yl 2 . 

Figure 01 shows the coupling constant c a evaluated for 
the harmonics in Table ||| We find that for all V (except 
very close to V c ) the NNN /-wave triplet is the dom- 
inant channel, while the other ansatze are repulsive or 
become repulsive upon increasing V. The performance 
of the NNN /-wave is understood by looking at the bot- 
tom panel of Fig. [3 The three marked points 120° from 
each other are nearly connected by x y/3 ordering 
wavevectors and the scattering among these points be- 
comes attractive upon approaching the CDW regime. As 
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Ansatz gk 


Label 


cos fc ■ ei + cos fc ■ e2 + cos fc ■ 


s 


sin k ■ e.% — sin k ■ e% 


Px 


cos fc ■ e.2 — cos fc ■ e$ 


dxy 


sin fc ■ ei — sin fc • e2 + sin fc ■ e3 


f 


sin fc ■ (ei + &2) — sin fc ■ (e2 + 63) + sin fc ■ (e3 — ei) 


NNN / 


cos fc • (2ei + e2) — cos fc ■ (ei + 2e2) + cos fc • (2e2 + e3) 
— cos fc • (e2 + 2e3) + cos fc ■ (2e3 — ei) — cos fc ■ (e3 — 2ei) 


i 



TABLE I: Triangular lattice harmonics evaluated in Fig. 2] 
ei, e2, e3 = e2 — e\ refer to unit triangular lattice vectors. 
The labels are to be taken as descriptive only. 



displayed in Fig.^ the NNN /-wave ansatz has its posi- 
tive lobes oriented precisely in these three directions and 
is able to utilize this attraction. 

Here, we note that generic next-nearest-neighbor 
ansatze utilize this nesting better than the nearest- 
neighbor ones, and we tried the corresponding NNN ver- 
sions of the s, p, and d-wave; the coupling constants are 
improved (not shown), but are still far from approach- 
ing the NNN /-wave ansatz. This is because each ansatz 
also has to pay the cost of repulsive scattering on large 
non-nested portions of the Fermi surface, and the NNN 
/-wave ansatz appears to be best here as well. For exam- 
ple, the coupling constant for the s-wave remains nega- 
tive and appears on the scale of Fig. 0] only very close to 
the critical point, despite the fact that it gains from all 
attractive V ee (0,9') < 0. 

In fact, we also solve the full eigenvalue problem spec- 
ified by Eq. H3U|) . and plot the first (maximal) and the 
second eigenvalue in the same figure. We find that the 
NNN /-wave coupling constant coincides with the max- 
imal eigenvalue in the entire range of V < V c (except 
maybe very close to V c ). The second eigenvalue is well 
separated from the first; the corresponding eigenvector 
over a large range of V does not have one of the sim- 
ple s, p, d or i character, but instead has eight lobes in 
momentum space, four of each sign. 

Finally, Fig. El shows the maximal eigenvalue as a func- 
tion of V for several dopings. In each case, the maximal 
eigenvalue corresponds to the NNN /-wave ansatz. The 
nesting displayed in the bottom panel of Fig.[3]that favors 
the NNN /-wave improves as we lower the doping, since 
the Fermi surface passes through the end points of the re- 
duced BZ when x = 0.194 (see also Fig. From Fig. [5] 
we see that the enhancement in the coupling constant 
is strongest and over the broadest range for x — 0.20. 
Observe also that the residual interaction from the no- 
double-occupancy constraint only (V — 0) also favors the 
discussed NNN /-wave channel on the triangular lattice. 

We conclude by discussing what these results mean for 
the scale of superconductivity. We see that because of 
the nesting the maximal c a can be relatively large com- 
pared to similar predictions in the non-nested cases (such 
as square lattice at finite doping, Ref. 0). However, if 
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FIG. 4: Coupling constants for the triangular lattice harmon- 
ics m Table □ (the s-wave ansatz is below the bottom of the 
plot). The doping is x — 1/3, and the critical V c = 2.11 sets 
the right-hand plot boundary. We also show the maximal and 
the second eigenvalue; the NNN /-wave coupling constant co- 
incides with the maximal eigenvalue for all V shown. 



0.20 




FIG. 5: Maximal eigenvalue for dopings x = 0.20, 0.30, and 
0.40. The critical V c from the quadratic fluctuation analysis 
is indicated by vertical line in each case, while the lst-order 
transition is indicated with the corresponding arrow near the 
bottom of the graph (cf. Fig. |5J. 



we use Eq. (I31f) literally, the obtained scale of T c is still 
very small: for example, if we take c a — 0.1, which is 
fairly large, then exp[— l/c a ] ~ 5 • 10~ 5 , which is tiny. 
Here we remark that one specific ingredient that can en- 
hance T c is missing in the present treatment, namely, the 
enhancement of the effective mass by the nearest neigh- 
bor repulsion while remaining in the Fermi liquid. More 
precisely, at the presented level of analysis, the effective 
hopping t c ff = xt is not renormalized by V, whereas it 
can become several times smaller when additional short- 
range correlations are included as discussed in Ref. y and 
in the following section. Since l/c a ~ t e s/V e s, this can 
have dramatic effect on the calculated T c . 
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IV. JASTROW-GUTZWILLER 
WAVEFUNCTION STUDY. DISCUSSION 

We now consider in some detail how the above results 
apply to the original spin-1/2 model. To this end, we 
have performed a systematic trial wavefunction study 
of the tV model. Specifically, we consider a family of 
Jastrow-Gutzwiller wavefunction o?! 1 ? 

VMnau . . . ) = e-Ei<i det [ii> a (Rj)] det fyM)] 

where {R} and {R'} denote the positions of spin-up and 
spin-down fermions respectively, and the wavefunction is 
nonzero only when the two sets do not overlap, which is 
the result of Gutzwiller projection. {ip a } refer to appro- 
priate single-particle states that are occupied in the "pre- 
projected" wavefunction. Each configuration of fermions 
is also weighted by a Jastrow factor defined via two- 
particle pseudopotentials u(r,— Tj), which puts additional 
correlations into the wavefunction. 

The Fermi liquid state is obtained by occupying ap- 
propriate plane wave states ipk, \k\ < hp. From our 
studies of the tV model, we conclude that already the 
nearest-neighbor Jastrow factor gives good control over 
local correlations and is sufficient for an accurate energet- 
ics distinction between the Fermi liquid and competing 
charge ordered state. 

As discussed at length in Ref. 0, lattice gas system 
described by the nearest-neighbor u(ri,rj) = W under- 
goes a transition to a V3 x \/?> state for large W . Our 
direct optimization with such single-parameter wavefunc- 
tion showed that this state is driven into the charge- 
ordered regime for sufficiently strong repulsion V and 
dopings in the range 0.27<x<0.50. We also mentioned 
that once this happens, we are no longer justified in us- 
ing plane waves for the preprojected orbitals. We should 
instead consider more general single-particle states; for 
example, we consider orbitals obtained by diagonalizing 
a "trial" Hamiltonian 

H C DW = - X« C L C > + ^AaCia , (33) 
ij i 

where we allow generic hopping amplitudes Xij an d site 
potentials <fi that follow the -\/3 x V3 ordering pattern. 
Specifically, we select the A sublattice and take xab = 
Xac = Xbc = f — k; k > makes the hopping more 
dice-lattice-like, while k < makes it more honeycomb- 
like. We also use the convention ipg = ipc = — <Pa/2 
and vary ipA- Observe that this structure of the trial 



Hamiltonian is suggested by the meanfield treatment of 
the Fermi liquid - CDW competition in Sec. [H] and in 
fact by any such meanfield. However, our wavefunctions 
have further local correlations built in by the Jastrow 
factors which allow the liquid to better accommodate 
the nearest-neighbor repulsion (also leading to significant 
renormalization of the bandwidth as discussed in Ref.0), 
while this is missing in the meanfield. 

We summarize our wavefunction studies for three dop- 
ings x = 0.24, 0.33, and 0.40. We report only the re- 
sult of the three-parameter W, k, ifA optimization. When 
k = 0, ifA = 0, we obtain the Fermi liquid state. As dis- 
cussed earlier, this restricted study should suffice for an 
accurate determination of the uniform liquid to CDW 
transition. 

For dopings x = 0.33 and 0.40, we find rather abrupt 
transitions to the \/3 x V3 state at V c (x — 0.33) ~ 3.5 
and V c (x = 0.40) « 4.0. On the other hand, for the 
doping x = 0.24, the system shows no ordering till very 
large V = 15 to 20 and maybe even higher. Compar- 
ing with Fig. [21 the behavior at this lowest doping is 
very different from the meanfield prediction. This is an 
extreme manifestation of the fact that the charges are 
able to effectively avoid each other and remain in the liq- 
uid state, which is captured by our Jastrow-Gutzwiller 
wavefunction, while the meanfield uniform state cannot 
accommodate this and pays large repulsion energy cost. 

Finally, having established the dominant local energet- 
ics, we also tried adding superconducting correlations on 
top of the renormalized Fermi liquid state. We tried sev- 
eral superconducting ansatz including extended s-wave, 
p + ip, d + id, /-wave, and also their NNN versions. How- 
ever, in our studies we were not able to detect any im- 
provement in the energetics upon adding superconduc- 
tivity. From this we conclude that if some such state 
appears at low energies, the condensation energy is still 
very small to be detected by direct numerical studies. In 
this situation, we are left to rely on approximate analyt- 
ical calculations such as discussed in the main body of 
this paper. 
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